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Abstract: Two analytical procedures (Crank’s method and Dincer’s method) for porous solid materials were reevaluated and used to 
determine moisture diffusion coefficients and moisture transfer coefficients for larch lumber subjected to drying. A diffusion-like equation 
was used to describe drying process data. The lumber was idealized in the modeling as infinite plates. The moisture transport process in¬ 
side the board was assumed to be one-dimensional. The macroscopic drying kinetics curves of larch timber at particular conditions were 
determined experimentally. Based on these data, calculation for both the moisture diffusion coefficients and moisture transfer coefficients 
by the Dincer’s analytical procedure were made. The dynamic moisture diffusion coefficients by the traditional Crank’s method were cal¬ 
culated. In general, diffusion coefficients calculated by the Dincer’s method were all higher than those by Crank’s method. These results 
could be due to the differences between two analytical methods and also different characteristics between solid moisture diffusion process 
and heat transfer process. Therefore the analysis and solution procedures of moisture diffusion differential equations need to be adapted in 
the future. With drying temperature’s increasing moisture diffusion coefficient (D) and moisture transfer coefficient (k) increases accord¬ 
ingly. Also the relationships between diffusion coefficients and temperature as well as material moisture contents were analyzed by using 
Arrhenius equation and bound water transport theory. 
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Introduction 

Kiln drying is a traditional operation method which has 
wide-range application in China wood industry; however, it is 
still difficult to analyze this industrial process in a satisfactory 
way. Wood drying is a process of simultaneous heat and mass 
transfer, and the mechanism of moisture transfer in wood is very 
complex due to material biological structure and variable drying 
medium parameters (temperature, humidity, velocity, etc). Many 
mathematical models can be used to simulate this drying process. 
The diffusion coefficient method provides a generalized way that 
works well for boards with moisture content below the fiber 
saturation point (FSP), (Simpson 1993; Miao et al. 2002; Cai 
2005). With diffusion approach, the moisture content gradients 
during drying process as well as the influences of process pa¬ 
rameters and material parameters on drying curve can be esti¬ 
mated, so that the goal of optimizing technical process and sav¬ 
ing energy can be achieved to some extent. 

The main objectives of this research are: (1) to measure drying 
kinetics curves of larch lumber(Zan'x gmelinii) under particular 
experiment conditions; (2) to detennine diffusion coefficients 
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and transfer coefficients with Dincer’s diffusion analysis; (3) to 
detennine diffusion coefficients as a function of moisture content 
by using Crank’s approach. Further objectives are to compare 
two kinds of information, illustrate the relationship between 
model predictions and variations in drying conditions and mate¬ 
rial characteristics. 

Experimental material and procedure 

Green larch lumbers were first put in a conditioning chamber 
with the temperature of (40±2) °C and a relative humidity of 
(96±2)%. After the conditioning parameters reached equilibrium, 
wood moisture contents in thickness direction of lumbers were 
kept below wood FSP. In other experiment chambers, the tem¬ 
perature was kept at 50, 60, 70, 85°C, separately, with relative 
humidity of (50±2)%, and air velocity was 0.5 m/s. Larch lumber 
was flat-sawn, free of visual defects, four to eight specimens 
were chosen for each test with specimens dimensions in 200 mm 
(width) Xl50 mm (length) x20 mm (thickness). Before experi¬ 
ments all specimens’ four edges were coated with two layers of 
epoxy to restrict moisture movement from four edge position, so 
that all diffusion activities were in wood transverse-radial direc¬ 
tion. 

Moisture diffusion test was conducted in experimental condi¬ 
tioning chambers, with control sensitivity temperature at ±1°C, 
and ±2% of relative humidity. A digital balance with sensitivity 
of 0.001 g and a weighing sensor with sensitivity of 0.1 g were 
used for the test. During experiments air temperature, specimen’s 
temperature, wet-bulb temperature, weight of all specimens 
could be recorded by a computer. After the specimens’ moisture 
distribution came to equilibrium with that of drying media con¬ 
ditions, the specimens were dried at (103±2) °C to constant 
weight, so that the oven-dry weight could be obtained. Before the 
analysis the following assumptions should be made: 
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1. The diffusion coefficient is constant in particular moisture 
content intervals, and the diffusion processes are one di¬ 
mensional, mainly in wood radial direction. 

2. The initial values of temperature and moisture content are 
unifonn within the specimens. 

3. The temperature and relative humidity outside surface of the 
specimens comes to equilibrium with that of the surround¬ 
ing immediately. 

4. The transport of heat, moisture, or gas proceeds in a manner 
which is symmetrical in space through out the bodies. 

Diffusion model solution procedure of two analysis 
methods 

Dincer’s method 


_ (M -M „ ) 


(5) 


The following dimensionless parameters are introduced for 
non-dimensional solution: 


Bi = 


k-X 

D 


( 6 ) 


Fo 


D ■ t 


(7) 



( 8 ) 


The transient moisture diffusion process in solid objects drying is 
similar to the heat transport process in these objects. The Fick 
diffusion law is exactly in the same form as the Fourier law of 
heat transfer, in which temperature and thermal diffusivity are 
replaced by moisture concentration and moisture diffusivity re¬ 
spectively. So, similar to the transient heat transfer analysis, we 
consider three situations for transient moisture diffusion. That is, 
the case the Biot number takes values as Bi<0.1, 0.1<Bi<100, 
and Bi>100 (Dincer 1998). For the case as Bi<0.1, there is a 
negligible internal resistance to the moisture diffusivity within 
the lumber, which is almost impossible in wood drying process. 
So it will not be considered here. Cases, such as 0.l<Bi<l00 and 
Bi>100, representing infinite internal - surface resistances and 
negligible surface resistances respectively, are the two common 
situations in wood drying practice. We therefore mainly consider 
these two cases in present analysis. 

The one-dimensional transient moisture diffusion equation in 
one-dimensional rectangular coordinate can be described by 
Fick’s second law, which has the following compact form: 


where, B ; is the Biot number for moisture loss, F 0 the Fourier 
number for moisture loss, and r the dimensionless distance. 

The simplified solution to the equation (1) is, 


E d = A - B 


(9) 


where 
A = exp 


( 0.2533 Bi) 
( 1.3 + Bi) 


B = exp (- /J 2 ■ Fo) 


( 10 ) 

( 11 ) 


The characteristic equation for root (u) is, 
fj — arctan(0.640443B/ + 0.380397) forO.KBidOO (12) 

// = fyj for Bi>100 (13) 


dM JL_{ D dM ) ( 1 ) 

dt dx ( dx J 

where, M is the moisture content of wood (%), t the time (h), D 
the diffusion coefficient (cnT-lT 1 ), and x is the space coordinate 
measured from the center of the board. 

Equation (1) is subject to the following initial and boundary 
conditions: 

^(x,0) = ^ = (M ; - - M e ) for 0.1<Bi<100, and Bi>100 (2) 

— (^(0,t)l = 0 for 0.1 <Bi<l00, and Bi>100 (3) 

dt J 

- £>^ -F<p(x ,t)) = k ■ tp(x ,t) for O.KBidOO, </>(X,t) = 0 
for Bi>100 (4) 

where, X=L(thickness of specimens) , (f) = (M — M ), Mj is the 

initial moisture content, and M e is the equilibrium moisture con¬ 
tent. 

The moisture content at a position of the solid wood is 
non-dimensionalized by the following equation: 


where, /(is root of the characteristic equation. 

From drying practice, we know that the dimensionless mois¬ 
ture content distribution can be regressed by an exponential 
function with the following form, 

E d = J ■ exp(- S ■ t) (14) 

By taking A=J, 

(jU 2 -Fo)=(S -t) (15) 

where S is the drying coefficient (1/h), E A the dimensionless 
moisture content for Dincer method and J is the lag factor. 

From the definition of the Fourier number, the moisture diffu¬ 
sion coefficient can be obtained in the following fonn: 

S ■ X 2 ) (16) 

B 1 ) 

Equation (16) can be used to determine the moisture diffusiv¬ 
ity values for lumber subjected to drying, with parameters (8i) 
and p being extracted from Eqs. (10) and (12) respectively. 

The equation determining the moisture transfer coefficients ( k ) 
can be written in the following forms: 
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k = 


f—1- 

" (1 - 3.94813 

In j)~ 

1a j 

(5.13257 ■ 

In J ) 


(17) 


The solution procedures for wood moisture diffusivity and 

moisture transfer coefficient models are as follows. 

1. The average moisture contents of lumber against time are 
determined experimentally and then non-dimensionalized 
using Eq. (5). 

2. The dimensionless moisture content are regressed against 
time in the exponential form of Eq. (14), the lag factor ( J) 
and drying coefficient ( S) are determined accordingly. 

3. The Biot number is detennined from Eq. (10) by using the 
lag factor. 

4. The root (ju) is calculated from Eq. (12). 

5. The moisture diffusion coefficient is determined by using 
Eq. (16). 

6. The moisture transfer coefficient is detennined by using Eq. 
(17). 

Crank’s method 


In this method, the governing model equation, initial and bound¬ 
ary conditions are similar to that of first method. By assuming 
the diffusion coefficient is independent of moisture content, the 
following equation for diffusion coefficients can be used (Siau 
1984): 

D = (E c f-4L c 2 (18 ) 

5 At 


where, D is the diffusion coefficients, L c is the half thickness of 
the board which is the diffusion direction, and t is the drying 
time. E c is the fractional changes in average moisture content at 
time (i), can be determined by following equation: 

(M - M ) (19) 

c (M e - M t ) 

According to Eq. (18), a plot of E c 2 versus ( t ) are linear with a 
slope of 5.1 D/L c 2 , so D can be calculated from the slope of a 
linear regression equation of the experimental data, D can be 
calculated as follows: 

D = 4I ^x Slope ( 2 °) 

5.1 


where, slope = ^ c ^ and the moisture content ( M) for the 
diffusion coefficient can be determined as 

M - M e +1 ( M; -Mj- 

Results and discussion 

Diffusion coefficients detennined by Dincer’s method (Dincer 
1995) 


matical regression analysis are shown in Fig. 1. For experimental 
conditions of 50% RH at 50°C, the lag factor and drying coeffi¬ 
cient are 1.06 and 0.0234 1/h, respectively, while the determina- 
tion coefficient ( R 2 ) is 0.891. As >7=1.06, this indicates that the 
Biot number is higher than 0.1, from Eq. (10) we know that B i 
equals 0.4049. So for this experiment the valid criterion is 
0.1<Bi<100. Then the root value (u) is deduced to be 0.569098. 
Having these data in hand, the moisture diffusion coefficient is 
found to be 0.072251 ctrr-h" 1 in Eq. (16) and the moisture trans¬ 
fer coefficient is determined to be 0.089219 cm-lf 1 by using Eq. 
(17). 

The experiment conditions’ parameters and the other experi¬ 
mental values are shown in Table 1 and 2, respectively. In Table 
2, the moisture diffusion coefficients have a increasing trend 
with condition temperature’s increasing, and the relationships 
between D and Temperature are similar more or less to that with 

( - ( 9,200 -70 M )3 

Arrhenius equation D = 0.07 -e RT (cal/mole), 

where M is moisture content and R is air universal constant and T 
is temperature in Kelvin. However, the detailed description of 
the relationship between D and T for larch wood need to make in 
future, which should be the objective of further research. 



Fig. 1 Measured and regressed dimensionless moisture content for 
larch board dried at 50"C and relative humidity of 50%. 


Table 1. The condition parameters of experiments 


Experiment no. 

Temperature( °C) 

Relative humidity (%) 

i 

50 

50 

2 

60 

50 

3 

70 

50 

4 

85 

50 

Notes: *—Dried from wood FSP; air velocity was 0.5m/s approximately. 

Table 2. Results of 

transfer coefficients 

moisture diffusion 

coefficients and moisture 

Moisture diffusion coeffi 
Experiment no. . 21 

cient (cm h ) 

Moisture transfer 
coefficient (cm-h 1 ) 

1 

0.072251 

0.089219 

2 

0.088445 

0.09695 

3 

0.093091 

0.052622 

4 

0.149985 

0.193067 


Diffusion coefficients determined by Crank’s method 


The dimensionless experimental moisture content distribution Moisture diffusion coefficients of different conditions are listed 

and the dimensionless moisture content variation from mathe- ' n Table 3. Because of the differences in calculation procedures, 
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the values of diffusion coefficients are all smaller than those of 
Dincer’s method, but the relationships between D and T in 
Crank’ method are similar to that of formers. With condition 
temperature’s increasing, the diffusion coefficients increase, too. 


Table 3. Average moisture diffusion coefficients calculated by 
Crank’s method 


Experiment no. 

Moisture contents 

(%) 

Diffusion coefficients 

(cm 2 -h J ) 

i 

18 

0.0106 

2 

16 

0.014 

3 

15 

0.0216 

4 

12 

0.026 



Moisture content(%) 


Fig. 2 Measured and regressed dynamic diffusion coefficients un¬ 
der condition with 70°C, 50% of RH 

The function relationship between diffusion coefficients (D) 
and moisture content ( M) is represented in Fig. 2. Below wood 
FSP, the moisture diffusion coefficients first decrease somewhat 
slowly and then decrease acutely in great speeds. This changing 
trend can be regressed very well by using polynomial function; 
the determination coefficient ( R 2 ) is 0.9931. The changing fonn 
under other experimental conditions is similar to the above one. 
The behavior of moisture diffLisivity with changing moisture 
content can be explained by theory of the bound water transport 
in wood. As we know, the activation energy (E b ) of wood bound 
water increases with moisture content decreasing. The approxi¬ 
mation expression of activation energy (E b ) is as follows. 

( 21 ) 


where, M is the moisture content. As a result, with moisture con¬ 
tent decreasing, more energy is needed to remove the bound 
water out of wood and difficulty for bound water movement 
increases accordingly. In macro level, the moisture diffusion 
coefficients decrease quickly with moisture content decreasing. 

Conclusions 

The moisture diffusion coefficients of larch wood lumber were 
calculated by using Dincer’s method and Crank’s method. Also 
the moisture transfer coefficients were detennined by using Din¬ 
cer’s method. The results of diffusion coefficients from the first 
approach were higher than that of the second one generally. 
These results may be due to the differences between two ana¬ 
lytical methods and also different characteristics between solid 
moisture diffusion process and heat transfer process. As a result, 
the solution procedures of relevant moisture diffusion differential 
equation need to be improved in the future. 

All the experimental values indicate that moisture diffusivity 
of larch wood lumber increases with increasing of drying tem¬ 
perature, partly testifies the validity of Arrhenius equation, how¬ 
ever the detailed study on relationship between D and Tempera¬ 
ture as well as Relative Humidity is future need. In order to gain 
thorough analysis of both calculated moisture transfer coeffi¬ 
cients and moisture diffusion coefficients which are determined 
by using the above two methods, dynamic moisture transfer and 
heat transfer tests with different specimen thickness are needed. 

Below FSP, as moisture contents decrease, Larch wood diffu¬ 
sion coefficients decrease quickly due to the increase of bound 
water activation. The relationship between D and M can be re¬ 
gressed very well by using polynomial function. 
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